弥散型核燃料元件热应力的计算模拟

针对弥散型燃料板的具体形状,选择了一种能够反映其实际变形特征的三维代表性体积单元作为研究对象;利用有限元法对其热力耦合行为进行了弹塑性分析,研究了热应力在燃料和基体中的分布特征.     

用迭代法模拟注塑成型的三维流动

Hele-Shaw模型是模拟注塑成型过程的常用模型,它的主要缺点是不能模拟一些重要的物理现象及引入实际制件中并不存在的“中面”概念,为了消除这些不足,本文开发了真三维的流动分析程序。建立了粘性、不可压缩的非牛顿流体流动的控制方程,为了避免同时求解耦合的压力场、速度场,本文引入拟稳定技术独立地求解这些方程,用迭代法求耦合方程的解。这种方法可以减少内存并提高数值方法的稳定性。算例表明数值结果与实验结果吻合较好,这种方法成功地模拟了注塑成型流动过程中的重要特征。     

On the bursting of linear polymer melts in inflation processes

Molten LLDPE and HDPE plates (thickness 2 mm) have been inflated into a circular cylinder (inner radius 31 mm) under isothermal conditions. Low deformation rates allow the plates to be inflated considerably into the cylinder, and at high inflation rates an early burst is observed.Axis-symmetric numerical simulation of the inflations have been performed, using a constitutive equation in the form of a separable memory integral where the strain dependence is described by the Linear Molecular Stress Function (L-MSF) model with dissipative convective constraint release. The material parameters in the constitutive model are obtained using liner viscoelastic (oscillatory shear) and uni-axial elongational measurements.The numerical simulations were performed for inflation of a flat plate and a perturbed plate, where a small circular cone was removed from the centre of the surface of the plate. This was done in order to investigate the stability of the inflations. It is shown that all of the inflations are hydrodynamically unstable, though the effect on the occurrence of the burst is limited. One exception is at slow inflation, where an unexpected burst may appear as a consequence of minute deviations from an ideal flat plate. All of the numerical calculations show quantitative agreement with the experiments for a wide range of experimental conditions. This strongly suggests that the initiation of the burst is a hydrodynamic phenomenon.The critical parameters in the inflation of molten linear polymers have been investigated using the Gel equation as a memory function (M(s)=Ans ^–(1+n)) and inflating the plate with a constant velocity for the top of the plate. The hydrodynamic burst in a linear polymer is mainly associated with the linear viscoelastic properties and only slightly with the non-linear strain dependence. Increased (linear) elasticity reduces the inflated volume, at the same inflation velocity, before the burst occurs. Furthermore, the critical parameter for the occurrence of the burst (whether or not the burst occurs) is related to the crossover point (G=G) in linear viscoelasticity.     

大型结构分析程序系统单元库更新

阐述了大型结构分析程序系统单元库更新的概念及其意义，报告了装入ＡＰ－１程序系统单元库的精化四节点及八节点平面等参元的情况，给出了更新前、后单元通过考题库测试的结果及与其它程序系统结果的比较。数值结果表明通过更新单元来提高程序系统的分析精度与计算效率从而增强其竞争力是可行的。     

中心集中力作用下圆薄板的固有频率—载荷特征关系

本文讨论了中心集中静载作用下圆薄板在非线性弯曲静平衡构形附近的微幅自由振动,其静平衡问题采用问题的精确解在此基础上,用伽辽金法获得了其最低固有频率——载荷特征关系.所得结果可供这类谐振弹性元件设计中参考.     

隐式Newmark法分析波动问题精度的探讨

本文研究了用平均加速度的Newmark积分格式求解一维平面弹性波动问题的精度。由于平均加速度的Newmark积分格式没有振幅衰减、不产生数值阻尼,因此采用此法分析波动问题会产生较强的高频振荡,尤其是加速度的计算值误差甚大。本文对防护工程中常见的平台型荷载和升压三角形荷载作用下的波动问题,应用有限元分析的实例表明:采用有数值阻尼的Newmark积分格式是改进有限元求解波动问题的一个有效方法,可以大大地提高应力和加速度求解精度。文中推荐的有数值阻尼Newmark积分格式的积分参数(=0.9,=0.49)可供工程中土与结构动力分析时参考使用。     

基于高阶位移模式的压电层合板动力有限元模型

建立了含压电片层合板的有限元动力学模型。以位于压电层上下表面处的电场强度和层间电压为未知量,给出了三次函数的电势分布模式,采用Reddy的高阶剪切理论描述板的位移场,假设板厚度方向的正应力为零给出了减缩的本构方程,采用有限元方法,基于Hamilton原理导出结构的动力学方程,然后用静态缩聚的方法压缩掉电场自由度和次要的位移自由度。最后用四边形矩形单元求解了一对称铺层和非对称铺层悬臂板的固有频率,并与ANSYS结果对比验证了本文模型的精确性。     

铁电薄膜在红外辐射作用下的响应

研究了红外探测器中红外焦平面列阵的象元———铁电薄膜微桥在红外辐射作用下的输出信号。用层合板壳作为微桥结构的力学模型,中间一层为压电材料,上下两层为金属材料电极。采用了力-热-电耦合的控制方程和变分原理,推导出了基于Mindlin假设的压电材料层合板有限元公式。以红外探测器在夜间从飞机上探测地面的坦克为例,用有限单元法模拟了铁电薄膜微桥在红外辐射作用下的力、热、电输出信号,并对结果作了分析比较。     

不确定非线性结构动力响应的区间分析方法

研究多自由度非线性不确定参数系统的动力响应问题. 以区间数学为基础，将不确定 性参数用区间进行定量化，借助一阶Taylor级数，给出了近似估计非线性振动系统动力响 应范围的区间分析方法. 从数学证明和数值算例两方面，将其与概率摄动有限元法进行了比 较，结果显示区间分析方法对不确定参数先验信息具有要求较少、精度较高的优点.     

Dynamic characteristic analysis of a 3-D semi-submerged body as a fluid-structure interaction system

An Arnoldi‘s method with new iteration pattern, which was designed for solving a large unsymmetric eigenvalue problem introduced by displacement-pressure FE (Finite Element) pattern of a fluid-structure interaction system, was adopted here to get the dynamic characteristics of the semi-submerged body. The new iteration pattern could be used efficiently to obtain the Arnoldi‘s vectors in the shift-frequency technique, which was used for the zero-frequency problem. Numerical example showed that the fluid-structure interaction is one of the important factors to the dynamic characteristics of large semi-submerged thin-walled structures.